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1. Introduction 

We study Wilson loop operators W^tf) in four-dimensional Euclidean SU(A'^) pure gauge the- 
ory, where 'rf is a. rectangular curve in M"*. Perimeter- and corner-divergences of W are eliminated 
by a (continuous) smearing procedure [1,2], where the associated smearing parameter s of dimen- 
sion length squared introduces an effective thickness for the curve 

In the infinite-A^ limit, the eigenvalue spectrum of the Wilson loop matrix exhibits a non- 
analyticity [1,3], separating a weakly-coupled short-distance regime from a qualitatively different 
strongly-coupled long-distance regime. At the transition point, the gap around -1 in the eigenvalue 
spectrum just closes. While smaller loops are insensitive to the compactness of SU(A'^), the full 
group is explored for larger loops (a key ingredient for confinement). The transition point, which 
depends on the shape of provides a natural scale for matching perturbation theory to the long- 
distance description provided by effective string theory [4]. Our goal is to determine, by numerical 
lattice gauge theory methods, how well loops of sizes barely on the strong-coupling side of the 
large-N transition can be described by effective string theory. 

Using the standard single-plaquette Wilson action, we have obtained Monte Carlo estimates 
for smeared rectangular- Wilson loops on a hypercubic lattice for various N's, couplings, volumes, 
and loop sizes from a database of 160 uncon^elated equilibrated gauge fields. All statistical eiTors 
quoted below are determined by jackknife with the elimination of one single gauge configuration 
from the set of 160 at a time. The range of couplings we use is 0.359 <b< 0.369, spaced by Ab = 
0.001 (the upper bound prevents spontaneous Z^{N) symmetry breaking [5] on all our volumes 
V > 12"*). Satisfactory statistical independence for our observables is obtained for gauge fields 
at neighboring b's being separated by 500 complete SU(2) updates combined with 500 complete 
over-relaxation passes. The N ^ °° limit is taken at fixed b = The set of A'^- values we use 
consists of N = 7, 11, 13, 19, 29. For the continuous smearing parameter, denoted by S on the 
lattice, we mainly use 0.2 < S < 0.4. The Wilson loops Wn on the lattice are defined by 

WN{Li,L2,b,S,V) = ^{TrY[Ui). (1.1) 

The product is over the links / in the order they appear when one goes once round a rectangle 
of sides Li_2- AH our fits are applied to 

WNiLuL2,b,S,V) = -logWNiLuL2,b,S,V). (1.2) 

When the loops are square, the two variables Li,L2 are replaced by one L with the understanding 
that Li,2 = L. 

In Sec. 2, we present our results for the large-A'^ string tension obtained exclusively from square 
loops. In Sec. 3, we extract a purely loop-shape dependent number from the data and compare it to 
the prediction made by effective string theory. For a more detailed presentation and discussion we 
refer to Ref. [6]. 

2. String tension from square loops 

We first want to determine \imN^c<, (limy^oo WN{L,b,S,V)) for square LxL loops. In principle, 
large-A^ reduction provides a shortcut for the infinite- volume limit. However, this requires tests and 
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fits since finite- volume effects depend onV,N, b and L. We use two different metliods to compute 
the limit: 

• Method 1) 

At fixed N, we compute wn on volumes that are sufficiently large for finite-volume effects 
to be negligible, then we determine WooiV = oo) (other arguments are omitted) by fitting 

Wn{V = oo) to 

= <») = »..(i' = <») + !!l<^ + 5^^. (2.1) 

We have evidence (strong for = 7 and = 1 1 , not that strong foiN= 19 and rather weak 
for = 29) that volumes V = 24^, IS'*, 14"*, 12"^ are sufficiently large for = 7, 1 1, 19, 29, 
respectively. This statement applies to the specific set of couplings and loop sizes we use. 

• Method 2) 

The second method makes use of large-A'^ reduction. At fixed V, we first take the limit N 
of W]\/{V) by fitting 

.MV)=.MV)^"-P^"-p. (2.2) 

So long as the center symmetry stays unbroken, there is no volume dependence in the infinite- 
A^ theory, i.e., Wcc(V) = w^{V = oo). We determine w„o(V = 12^) from A^ = 11, 13, 19, 29 
[method 2a)] and Woo{V = 14^) from A^ = 7, 11, 13, 19 [method 2b)]. 

We obtain reasonable values of X^/^dof for the fits^ and good agreement (compatible with the 
statistical accuracy of about 0.1%) between the three results for limN v^ooWN{V). Truncating the 
expansions (2.1) and (2.2) at 0'{N^'^) would result in very large X^/^dof, so aj cannot be set to 
zero. Including the N = 29 result in the fit (2.2) for V = 12^ is crucial for 2a) to agree with 2b) 
and 1) for large loops and large b. Including N = 7 in the V = 12^ fit would require an additional 
1 /N^ correction in (2.2). When the lattice size V 4 is getting close to the critical lattice size Lc{b) at 
which the center symmetry brakes, there is no useful information to be gained about the N,'V = °° 
limit from numbers obtained at low values of A'^. Since the required computation time scales as 
N^V , 2a) is about 1.75 times more expensive than 2b) and 1) is about 2.5 times more expensive 
than 2b). However, it is hardly possible to conclude from 2b) or 2a) alone that the estimates for 
Woo(V) are reliable. We became confident that we have correctly determined Xmvi^ v ^•^oW^iy) only 
after having obtained agreeing results from 1), 2a) and 2b). 

For fixed b and S, we use the shorthand notation Woo(L) = \imNy^ooWj^{L,b,S,'V) and expect 

w^{L) + ^\ogL^ = ci+C2L + oL^ + ff(^-^^ . (2.3) 

The log term comes from the determinant of small fluctuations around the minimal area configura- 
tion in the effective string description. We shall return to it in Sec. 3. For now, its presence is just 
assumed because including it gives good fits while excluding it gives bad fits. 



' 2a) is an exception where we have X^/^dof up to 6 for < 0.361 . This probably reflects the impact of the infinite-A' 
bulk transition at b = 0.360. 
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Neglecting corrections of order we fit 

1 (.4^ + 1) - »44 + ^ log (i + i) ) = (i. + + f + ^ (^) (2.4) 

to a straight line as a function of L + ^ to determine the string tension a and the coefficient of the 
perimeter term C2 (see Fig. 1 for some examples). Most 5x5 loops fall into the neighborhood of 
the large-A'^ phase transition in the eigenvalue spectrum of the Wilson loop matrix for the b and 5 
values we work with. Physically smaller loops will have a single-eigenvalue distribution which has 
a gap around - 1 . Therefore, we only use loop sizes in the range 6 < L < 9 to determine the string 
tension. 

Aw/2 




Figure 1: Plots of ^ = 5 (wc»(L+ 1) - Woo(L) + jlog (l + i)) obtained with method 1) as a function of 
L + 5 at 5 = 0.4 and /7 = 0.36 (red), b = 0.362 (black), b = 0.365 (blue) and b = 0.368 (green). Error bars 
are not visible in the plot. The straight lines show linear fits through the corresponding data points. Only 
points 6 < L + J < 9 are used in the fits. 



The results obtained for a using the different methods for computing Woo agree with each other 
within statistical errors, which are smallest for method 1). See Fig. 2 for a plot. As expected, a 
does not depend on the smearing level S within the en^ors (which increase with decreasing S). 

A scale length in lattice units denoted by ^c{b) is used to carry out extrapolations to the 
continuum limit. It is defined at A'^ = 00 using a three-loop calculation of the /3 -function for 
the lattice couphng. The coefficients are written as jSo = ^ = jSi = ^ = ji^^pji, and 
P2 = lim^f^oc ^ ~ —3.12211 X 10^^. Integrating the RG flow, we define: 

Above we have replaced the gauge coupling b by b]{b) = limAry^oo^W/v(L = l,b,S = 0,V), a 
substitution known as tadpole improvement. The definition of ^c{b) is taken to match with [7]. We 
only added a numerical prefactor to make ^c{b) ~ Lc{b), where (Z?) is given in [5]^. 

^This approximation is good to 10-15% in our range of couplings and would become exact at b = °°. 
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We separately carry out two two-parameter fits of the relation between the string tension (y{b) 
and ^c{b). The two pairs of parameters are denoted by do, dy and /o, f\: 



a{b) 



do 



+ 



d\ 



and 



1 



■ f^'a{b)+ha{bf 



(2.6) 



We use ranges 0.359 <b< 0.369 (range A) and 0.362 <b< 0.367 (range B). We also use the 
limited b range (B) since we have observed increasing ;^^'s for the infinite -A^,y extrapolations 
using method 2a) for b < 0.361, as mentioned above. Another reason is that finite-volume effects 
increase with increasing b. This reason only applies to method 1). The difference between the two 
fits is a simple indicator of systematic errors induced by the truncation of the perturbative series. 
We find that these particular systematic deviations are of the same order as the statistical errors. 

Our result for the infinite-A'^ continuum string tension is given by lim^^oo G{b)^^{b) = 1. 6{l) {3). 
The first error is statistical and the second systematic. The systematic error is more of a guess than 
a well founded estimate. In terms of A^, this translates to a/A^ = 3.4(2) (6). 
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Figure 2: Plots of CJ as a function of ^: method 1) in black, method 2a) in green, method 2b) in blue, 
together with corresponding fit functions (the fits are obtained using 0.359 <b< 0.369). 



Our central number is 2-3 standard deviations smaller than that of AUton et al. [7] (a/A^^2 = 
3.95(3) (64) a.tN = oo). Their numbers were extracted from Polyakov loops which are substantially 
longer than one side of our square loops. The systematic errors are dominated by the continuum 
extrapolation and their relative size is roughly the same for us. 

A previous estimate for the string tension at infinite A'^ extracted from rectangular Wilson 
loops has been given in [8]. Expressed in terms of our variables it is G^^\b,=o.i^2 = 2.2(3). While 
writing up our paper [6] a new study [9] appeared which also deals with rectangular Wilson loops. 
These authors obtain a/A-^^i = 3.63(3) (statistical error) at = oo if they apply the continuum 
extrapolation method of [7]. This number is fully consistent with ours and has very small eiTors by 
comparison. 

There seems to be a disagreement at the statistical level between [7] and our result which 
agrees with that of [9]. The result of [8] seems to side with that of [7], but has too large errors to 
be sure. The systematic en^ors ai^e too large to claim evidence for a difference between the string 
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tension extracted from Wilson loops and that extracted from Polyakov loop correlators, wliicli 
would be very difficult to accept at the theoretical level. 

3. Shape dependence 

We now turn to a study of the shape dependence of the scale-independent term in wn and 
compare it with the effective-string prediction. For rectangular Li x Lj loops it is convenient to 
introduce the modular invariant shape parameter 

The accuracy we now need does not permit taking the N ^ °° limit. We therefore restrict our 
attention to the N = 1 ,\ \ data. We shall see that the numbers we compute are identical within 
errors for N = 1 and N = \ \, indicating that it is unlikely that they will change in a substantial 
manner in the N = °° limit. 

At fixed b, S, V, and fixed finite N, we expect (arguments b, S, V are omitted) 

wn{Li,L2) + I logLiLz = ci,;v(C) +C2./v ^' + ONUL2 + — ]---] . (3.2) 

After having determined the lattice string tension and the coefficient of the perimeter term 
C2,N from square L x L loops (using the method described in Sec. 2) at fixed A'^, b, S, V, we fit 
wn{L,L) + ^ logL^ — OnL^ — C2.^L to a constant, ci a'CC = 2) (using loop sizes 6 < L < 10). Next, 
we analyze the results obtained for a sequence of rectangular loops at the same b, S, V, N with 
L2 = 2Li, i.e., C = i fixed. Using the results for On and C2,n obtained from square loops, we then 
determine ci,a?(^ = §) by fitting wn{L,2L) + ^log (2L^) — OnIL^ — C2,n^L to a constant (using 
4 < L < 7). Figure 3 shows a plot of ci Af=ii (2.5) — ci.a?=ii (2) as a function of b. Within statistical 
errors, our results for ci,a?(2.5) — ci^Ar(2) do not depend on Zj, S, orN. The effective-string prediction 
for ci (2.5) — ci(2) is 

ilog(5(^).-0.08«4. (33, 

where ri{x) is the eta-function. We find that the effective-string prediction is smaller than the 
observed values by a factor of about 1.5 to 1.7. 

We use sequences of L x L it 1 loops to cross check our results for the string tension and the 
shape dependence of c\^n. Furthermore, we simultaneously fit L x L, L x 2L, and L x Lit 1 loops 
to the functional form (3.2), where we expand ci jv around ^ = 2 and allow the coefficient of the 
logLiL2 term to become a fit parameter. This analysis confirms the expected value of 1/4 for the 
coefficient of the log term as well as the results for Gn and ci,n presented above. 

This means that one prediction coming from the determinant of Gaussian surface fluctuations 
in the effective string description works close to the large-A'^ transition in the eigenvalues and the 
other does not. However, these two predictions are somewhat different even within effective string 
theory (see [6] for details). 

Our result for the shape dependence of the scale-independent term is in agreement with [9] 
who independently report a deviation from effective string theory. 
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Figure 3: Plot of ci,a,(2.5) -ci.a,(2) for = 11 (on V = 18"^) as a function of for S = 0.2 (red), S = 0.28 
(green) S = 0.4 (blue), and S = 0.52 (black). 

A detailed discussion of possible explanations (from perturbation theory and by higher-order 
terms in the string expansion) is presented in [6]. Our data is not conclusive enough to settle this 
issue. We think that the shape dependence of planar Wilson loops presents an interesting case for 
testing the limitations of the effective string approach which deserves further study in the future. 
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